CMB Lensing — Power Without Bias 
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We propose a novel bias-free method for reconstructing the power spectrum of the weak lensing 
deflection field from cosmic microwave background (CMB) observations. The proposed method is in 
contrast to the standard method of CMB lensing reconstruction where a reconstruction bias needs 
to be subtracted to estimate the lensing power spectrum. This bias depends very sensitively on the 
modeling of the signal and noise properties of the survey, and a misestimate can lead to significantly 
inaccurate results. Our method obviates this bias and hence the need to characterize the detailed 
noise properties of the CMB experiment. We illustrate our method with simulated lensed CMB 
maps with realistic noise distributions. This bias-free method can also be extended to create much 
more reliable estimators for other four-point functions in cosmology, such as those appearing in 
primordial non-Gaussianity estimators. 



When the universe was ^380,000 years old, the pho- 
tons of the cosmic microwave background (CMB) decou- 
pled from the primordial photon-baryon fluid. Travel- 
ing through the newly transparent universe, the photons 
were deflected many times by the gravitational influence 
of the large-scale structure potentials through which they 
passed, an effect known as the gravitational lensing of the 
CMB, or CMB lensing (for a review, see [H). The net 
effect of lensing is an effective deflection of each CMB 
photon, leading to a remapping of the points on the 
CMB sky. The deflection field depends on a distance- 
weighted projection of density perturbations along the 
line of sight. The power spectrum of the deflection field 
is therefore sensitive to both geometry and the growth 
of structure over a broad redshift range (z ~ 0.5 — 5). 
As such, knowledge of the convergence field can provide 
strong constraints on parameters that affect geometry 
or the growth at later times, such as the sum of neu- 
trino masses and parametrizations of non-standard dark 
energy behavior [2j. These constraints are complemen- 
tary to those obtained directly from the primordial CMB 
anisotropics. In the very near future, ongoing and up- 
coming CMB experiments, such as Planck, ACT, SPT, 
PolarBear, ACTPol and SPTPol will produce datasets 
with sufficient resolution and sensitivity to begin the de- 
termination of the deflection field and ultimately realize 
the cosmological potential of CMB lensing science. Ro- 
bust algorithms that are insensitive to the details of the 
noise properties of the survey will be essential for accu- 
rate determination of the deflection power spectrum. 

In this Letter, we propose a novel and bias-free tech- 
nique for the measurement of the lensing deflection power 
spectrum. This is in contrast to the standard optimal 
quadratic estimator (OQE) method where a bias term, 
comparable to and often much larger than the signal, has 
to be subtracted (for details, see 0,3)- This bias, which 
is a temperature-field four point function that depends 
on noise and foregrounds, must typically be computed 



or simulated to an accuracy of a few percent in order to 
get reliable estimation of the signal. Because of our lim- 
ited knowledge of the temperature and polarization fore- 
grounds, and the typically complicated noise properties 
of CMB experiments, modeling this bias term sufficiently 
accurately for a robust detection of lensing may not be 
possible. 

To discuss how this bias appears and can be avoided, 
we will review some lensing theory here. The lensed and 
unlensed temperature fields (unlensed quantities will be 
denoted by a tilde) are related by T(n) = T(n + a(n)), 
where a(n) is the deflection field and k = • a is the 
convergence field. In the flat sky approximation, the 
temperature field can be expanded as a Taylor series in 
Fourier space @, @] : 

r rl 2 f' £'■(£- £') 
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Thus, gravitational lensing introduces correlations be- 
tween the formerly independent modes of the temper- 
ature field, which can be used to construct a quadratic 
estimator for k: 



k(L) = N k (Jj) 
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where N K is a normalization that ensures that the esti- 
mator is unbiased (i.e., (£} CMB = k) and Fw and Fa 
are filters that can be tuned to minimize its variance |3| . 
Here we will use notation from Q and denote ensem- 
ble averages over CMB realizations with the large scale 
structure (LSS) fixed as (...) CMB , and averages over LSS 
realizations as (...) L gg, and use (...) to denote ensemble 
averages over both, i.e., ((■•■)cmb)lsS" 

While equation ([2]) provides an unbiased estimate of 
the convergence field, the naive lensing power spectrum 
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estimator k*(L)k(L) is highly biased: 
(k*(L)k(V)) = N K *(L)N K (L') 

r g(£,L)g(£',V 
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where L) = L • (L - £)F W (£)F G (\L - In going 
from the first equality above to the second, we have ne- 
glected a few terms that appear involving integrals over 
the convergence power spectrum. These higher order 
terms are computed in [J] and are subdominant com- 
pared to the Gaussian four point term above. Applying 
Wick's theorem contractions to the four-point term, the 
above equation can be reduced to: 

(/s*(L)ft(L')) = (2tt) 2 (5(L - L')(C£ + c£' Gauss ), 

where 
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with /(.£, £, L) ~ g(£,h) g(£' This expression shows 
that the Gaussian bias (q£< Gauss depends on the map 
power spectrum estimate which is usually a sum of the 
CMB temperature power spectrum, foregrounds, and 
noise. Since the Gaussian bias term is typically more 
than an order of magnitude larger than the intrinsic lens- 
ing signal C^, the standard approach to estimating lens- 
ing requires a detailed understanding of each of these 
contributions. 

The goal of this paper is to eliminate the need to com- 
pute this Gaussian bias term at the percent level by elim- 
inating the bias altogether. From (|3|) and ([5|), one can 
see that the brute force way to achieve this would be to 
perform the double two-dimensional integral explicitly in 
(|3]) with the following conditions imposed on the function 
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if £' = £ or £' = L 



f(£, £' ', L) otherwise. 
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However, this method is computationally expensive and 
not efficient. A conceptually simpler and much more ef- 
ficient method of eliminating the bias would be to parti- 
tion the Fourier space into non-overlapping annuli, and 
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FIG. 1: Graphical expansion of the naive estimator after split- 
ting up the Fourier space into an inner and an outer annulus. 
We use the linearity of both operators in this expansion. The 
terms with Gaussian bias are shown enclosed by boxes. The 
underlined terms (identical by symmetry) are implemented 
in the simulations described in this paper to illustrate the 
method. 



cross correlate the k(L)'s reconstructed from tempera- 
ture modes in two disjoint annuli. To formulate this fully, 
let us introduce some compact notation. Let us denote 
the operation of reconstructing n from two temperature 
maps, equation ©, (T (£> T) where all quantities 

are understood to be in Fourier {£) space. Now, the naive 
estimator for the convergence power spectrum ([3]) can be 
written as k'k = (T®T) x (T<g>T). Now consider break- 
ing up T(£) into two Fourier space maps, one which has 
non-zero elements only with an annulus £q < I < £i (the 
"in-annulus" ) and another which has non-zero values in 
£i < I < £ 2 (the "out-annulus" ), where £ < £\ < £ 2 - 
Writing T{£) = T in (£) + T out (£), the naive power spec- 
trum estimator can be written out as, 
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(T T) x (T <g> T) 

= ((T in +T out )® (T in + T out )) 

x((T m + T out )®(T ln + T out )). 



(7) 



We have expanded this out graphically in Fig. [T] In this 
figure, we represent the in-annulus by a filled circle and 
the out-annulus by an empty circle. The expansion leads 
to 16 terms, some of which are identical due to symme- 
try. Note that the Gaussian bias associated with each 
term evaluates to a sum of two terms from two possible 
Wick's theorem pairings (one member of the pair being 
taken from either side of the "x" sign). Those terms 
(e.g. (T in ® Ti n ) x (Tout <8 T out ), underlined in Fig. □ ) 
where either pairing leads to at least one product of an 
out-annulus with an in-annulus will have contributions 
to the bias which evaluate out to zero. We find 10 such 
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FIG. 2: Left: Convergence power spectrum reconstructed with the proposed Gaussian bias-free method from four 5° x 15° 
patches with simulated CMB signal with 2 pK-arcmin white noise. The blue (filled) circles show the mean of 120 Monte Carlo 
realizations with lensed CMB, while the green (empty) circles show the same for unlensed CMB maps. The error bars are 
estimated from the scatter between Monte Carlo runs and are representative of the uncertainty expected in one realization; the 
lensed errors are higher than the null errors due to the presence of a sample variance component. The red continuous curve 
is the input theory for the convergence field power spectrum. Right: Same as left, but for non-white and anisotropic noise 
simulations seeded by noise in ACT maps, reduced in amplitude by a factor of 3. 



terms in Fig. [TJ the remaining 6 terms' biases have non- 
zero expectation value and hence contribute to the total 
bias (these terms are shown enclosed by boxes). There- 
fore, one can construct a new estimator for the conver- 
gence power spectrum by applying these annular filters 
and optimally combining the 10 terms (some of these are 
identical) that are bias-free. Of course, eliminating bias 
comes at the cost of reducing the signal-to-noise because 
we throw out a fraction of terms with information. In 
principle, higher signal-to-noise can be achieved by fur- 
ther subdividing each annulus into an inner and an outer 
part, and iterating the method, thereby reducing the ra- 
tio of the number of biased to bias-free terms. 

There are a few details that need to be taken into 
account when applying this method to a real experi- 
ment. For a partial sky map, nearby Fourier modes will 
be coupled by the power spectrum of the data-window 
with some characteristic width A£. There can be addi- 
tional effects such as coupling of nearby modes induced 
by anisotropic noise. In general, if the effective width A£ 
of such correlations is known, the above method needs 
to modified by separating the two annuli by some small 
multiple of At. This will ensure that our annular method 
for eliminating the Gaussian bias works despite these cor- 
rellations due to systematics. 

Due to the annuli being separated by A£, all the terms 
involving a convergence map obtained from an innner 
annulus as well as an outer annulus are undefined for 
£ < A£ (as can be deduced from equation^. Hence, 
for a simple bias-free estimator we used only the terms 



underlined in Fig. [TJ those with one "in-in" convergence 
map crossed with one "out-out" convergence map, so 
that the new bias free estimator is: 

^bias-free = ^ g ^ ^ ^ ^ ^ (g) 

There is one further subtlety to the reconstruction of 
the convergence power spectrum. Window function cor- 
relations due to finite maps and anisotropic noise not only 
affect the Gaussian bias, but also appear more directly 
as a spurious lensing signal in the reconstructed conver- 
gence maps. This spurious convergence, (k), needs to 
be simulated and subtracted off from the reconstructed 
k map. To determine whether the reconstructed con- 
vergence power is sensitive to the accuracy of the sim- 
ulation of (k), we must determine the magnitude of 
C/T ^ out cmt (where one must distinguish between 
the two (k) fields because they are simulated with dif- 
ferent annular filters). From our Monte Carlo simula- 
tions, which we describe below, we find the following: 
for our bias-free lensing estimators, ,n, * n ^°" t °" t ap- 
pears consistent with null and is typically two orders of 
magnitude smaller than the true convergence power CI 
and three orders of magnitude smaller than the previ- 
ously discussed reconstruction bias. We also verified that 
the results of our simulations plotted below are the same 
whether or not we subtract (k) from k. The accuracy of 
the simulation of the spurious signal (k) thus seems to 
only have a negligible influence on how well the lensing 
power spectrum is reconstructed. 
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We illustrate our method for bias-free lcnsing power 
reconstruction by performing Monte Carlo simulations, 
loosely modeled on the observations made by the Ata- 
cama Cosmology Telescope (ACT). We choose our sur- 
vey geometry to be an oblong 5° x 60° stripe, divided 
into four adjacent patches of 5° x 15° each. We simu- 
late convergence maps on these patches from an input 
power spectrum, and generate deflection fields from the 
convergence maps. We also generate Gaussian random 
realizations of unlensed CMB maps from an input power 
spectrum, which we subsequently lens using the simu- 
lated deflection field. Then we smooth the maps with 
the ACT beam (1.4 arcmin full- width-half-maximum), 
and add noise. We perform two variations on the noise. 
In the first version, we add white noise at the level of 
2 u-K-arcmin. In the second version, we simulate noise 
seeded by the noise power spectrum realized in the maps 
from ACT. These simulations capture the non-white and 
anisotropic aspects of the noise in ACT maps (for the 
detailed procedure see Q). To ensure that the recon- 
structed convergence power spectrum is not too noisy 
(so that the Monte Carlo simulations rapidly converge), 
we reduce the amplitude of the simulated noise by a fac- 
tor of three over what is observed in ACT maps. These 
maps are roughly 10 uK-arcmin in noise. For each type 
of noise, we simulate 120 realizations of the full map by 
randomizing the CMB, the convergence and the noise. 
We then apply the bias-free convergence power spectrum 
estimator method to each random realization of the noisy 
maps. We apply this both to lensed maps with noise and 
as a null test also to unlensed maps with noise. We de- 
fine our annular filters such that the inner annulus is 
I = (500 - 1500), the outer annulus is I = (1900 - 3300), 
so that A£ = 400. 

The results are shown in Fig. [2j This figure shows 
that in both cases our method is able to extract the con- 
vergence power spectrum without bias (note that a small 
amount of bias from higher order terms discussed in Kes- 
den et al. Q is present in the reconstruction) . The figure 
also shows that with noisy but unlensed CMB maps, the 
measurements are consistent with a null signal. 

Finally, we should point out some caveats. Here we 
have assumed CMB as the only signal. In reality, emis- 
sion from point sources and the thermal and kinetic 
Sunyacv-Zcldovich effects contribute to mm-wave maps. 
Also, noise correlations in real experiments can be more 
complicated than what is simulated here. There are also 



other subdominant sources of reconstruction bias, such as 
those discussed in Q and Q . More work will be needed 
to characterize these foregrounds and biases in context 
of our new method. These will be discussed in a future, 
more detailed work. 

It should be noted that our approach for bias-free lens- 
ing convergence reconstruction can be easily extended to 
estimating other four-point-functions in cosmology, such 
as the estimators of primordial non-Gaussianity. 
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